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Abstract. The problem of finding independent components of an indexed object (e.g., a 
tensor) with arbitrary number of indices and arbitrary linear symmetries is discussed. It 
is proved that the number of independent components f(k) is a polynomial of degree not 
greater than the number of indices n, k being the dimension of the space. Several algorithms 
to compute f(k) for arbitrary k are described and discussed. It is shown that in the worst 
case finding f(k) for arbitrary k requires solving at most P(n) systems of linear equations 
with at most (n!) 2 equations for at most of n! unknowns, P(n) being the number of partitions 
of n. As a by-product, an efficient algorithm to parametrize all components of the object 
through its independent components is found and implemented in Mathematica. 
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Q\ ' 1 Motivation 

In this paper we attempt to find an approach to the problem formulated, e.g., in [3]. The problem 
consists in calculating the number of independent components of an indexed object A y with 
symmetries (for example, A 1 ^--- = A i'"). The importance of this problem for computer algebra 
systems for indicial tensor computations can be clarified by the following example given in [3] . In 

(\ | many packages for indicial tensorial calculations one can define an object A u with 3 indices and 

Q H . the following two symmetry properties: 

feb; A ijk = A jik (i) 

and 
k^ ; A i: > k = -A ikj . (2) 

One can easily check that these two symmetry properties imply that all components of A vk is zero. 
However, none of the packages available to the author can automatically recognize this fact. In many 
cases, however, if the software could recognize such and similar situations, it could simplify many 
kinds of calculations. On the other hand, it is often quite important to find not only the number of 
independent components, but also the independent components themselves (for example, to store 
the components of an object in the best possible way). In this paper we will give efficient algorithms 
for both finding and counting the independent components of an object with symmetries. 

2 The Problem and Some Trivial Properties 

Here the following problem with linear symmetries is considered. Let A lll2 '" in be an object with 

n indices. Each index ij can take values from ij G {1, 2, . . . , fc}, where k > 1 is the dimensionality 

of the space where the object is defined. The object has s symmetries defined by equations of the 

form 

pi 

J2ajiA*' l =bi, I = 1,2,... ,a (3) 

where a,ji and bi are numbers characterizing the l-th symmetry, and %ji is a permutation of the n 
indices i\i2 ■■■in, Pi is the number of terms in the l-th. symmetry. Clearly, the maximal number 
of the terms in one symmetry and the maximal number of linearly independent symmetries are 
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constrained as (otherwise the individual terms and the whole symmetries linearly depend on the 
other terms and symmetries, respectively) 

Pi < n\ , (4) 

s < n! . (5) 

In principle, the number of symmetries s can be larger than n\, but in this case either the additional 
symmetries are linearly dependent on the first n\ symmetries or at least one of the additional 
symmetries are incompatible with the first n\ ones and such an object does not exist. A component 
of A lll2 '" ln is called independent if it is not constrained by symmetries (3) to be a number and 
if it cannot be calculated as a function of other independent components. The set of independent 
components is obviously not unique (if a symmetry requires A 21 = A 12 any of these two components 
can be considered as independent). The size of this set is, however, independent of this non- 
uniqueness, and it is sufficient to find any version of it. The problem consists in computing of the 
number / of independent components of A ll%2 '" %n as a function of k: f(k). Since the total number 
of components of A lll2 '~ ln for a fixed k is k n , one concludes 

< f(k) < k n . (6) 

3 Objects with One and Two Indices 

In order to understand the problem more clearly, let us consider first all possible objects with one 
and two indices. The case of n = 1 is quite trivial. An object A 1 has at most k components. If the 
object has no symmetries of the form (3), the number of independent components f(k) = k. The 
only possible symmetry reads 

aA i = b, (7) 

where a and b cannot vanish simultaneously (otherwise the object has again no symmetry). Clearly, 
if a = and b =/= the symmetry cannot be satisfied. Such symmetries will be called incompatible. 
In case of incompatible symmetries we will symbolically write f(k) = 0. If a =/= then A 1 = b/a, 
i.e. all components for any k are constrained to be number b/a. This can be summarized as follows 



f(k) = 



k, no symmetries (a = b = 0), 

0, a^O, (8) 

0, a = & b j- 0. 



The case of n = 2 is a bit less trivial. For an object A*i with no symmetries f(k) = k 2 . Each 
of the two possible symmetries is of the form 

aiA ij + a 2 A ji = b. (9) 

Let us consider first only one symmetry (9). For C\ — k "diagonal" components with j = i one 
can write this symmetry as (no implicit summation is assumed here) 

(oi + a 2 )A u = b, * = l,2,...,fc. (10) 

This symmetry can be analyzed in the same way as we did for (7) above. The k "diagonal" 

components are constrained independently of each other (but clearly in the same way). Therefore, 

it is sufficient to analyze (10) for one fixed i and then multiply the result (either or 1 independent 

components) by k. Among the other k(k — 1) "off-diagonal" components with j =/= i there exist 

a! 
C 2 pairs of components, C h a = — — being the binomial coefficient. It is only the components 

within each of these pairs which could be potentially constrained by (9). Indeed, for any fixed j and 
i (j 7^ i) only two combinations of indices, ij and ji, are connected to each other by a permutation 
and appear in (9). Therefore, it is sufficient to consider any of these pairs and multiply the results 
(one can have 0, 1 or 2 independent components within each pair) by C\ = k{k — l)/2. Therefore, 
the number of independent components of A tJ for an arbitrary k can always be calculated as 
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f(k) — f 2 k (k — l)/2 + /i k with two integers < f 2 < 2 and < /i < 1 depending on the 
symmetries of the object. Let i — 1 and j = 2. Then (9) give two linear equations constraining A 12 
and A 21 : 

ai A 12 + a 2 A 21 =b, 

ai A 21 +a 2 A 12 =b. (11) 

Analyzing (10) and (11) one gets the number of independent components of an object with n = 2 
and one symmetry (9): 

k 2 , no symmetries (ai — a 2 = b = 0), 

i fc (fc + 1), ai = -a 2 + & 6 = 0, 

f fc(fc-l), ai=a 2 ^0, (12) 

0, |ai|#|a 2 |, 

0, ai + a 2 = 0&6^0. 



/(^ 



Clearly, for ai = — a 2 ^ & fo = the symmetry (9) can be written as A*i — A> 1 = 0. This means 
that the "diagonal" components arc not constrained at all, and the "off-diagonal" components are 
pairwisc equal. We have thus a symmetric "matrix" . For a\ = a 2 ^ one has A*i + A 11 = b = 
b/a\. This means that the "diagonal" components are all equal to 6/2 and the "above-diagonal" 
components can be computed from the "below-diagonal" ones as A^ = —A> 1 + b. For b = we 
have a skew-symmetric "matrix" here. In case |ai| ^ |a 2 | all components are equal to the same 
number: A % i = b/(ai + a 2 ). 

It is easy to see that for objects A % i with two symmetries 

an^ ij +a 2 iA^=6i, 

a 12 A^ + a 22 A ji =b2 (13) 

f(k) also takes the same 5 possible values as in (12). The ideas found in these examples will allows 
us to formulate several important theorems and algorithms below. 

4 The Number of the Independent Components 

Let us first leave aside the question whether the symmetries of the object are compatible and 
consider that f(k) ^ for any k. This will be further discussed in Section 5. The results of the 
previous Section suggest to formulate the following 

Theorem 1. For an object A llt2 •••*" with arbitrary n and arbitrary compatible linear symmetries 
the number of independent components f(k) is a polynomial of the dimensionality k. The degree 
of this polynomial does not exceed n. 

Proof. From (6) it is obvious that if f(k) is a polynomial its degree cannot exceed n. 

Let us denote the set of all possible values of each of the indices i a as Nk — {1, 2, . . . , k}. Its 
cardinality is \Nk\ = k. Let us also denote the set of all independent components of an object 
A ili2 - in with all i a e N k as E{N k ). The cardinality of E(N k ) is clearly f(k). It is clear that we 
can change N k to any set S with the same cardinality k and number of independent components 
of 4»i»2...»„ with ia e 5- wi n be again f(k). That is, for any S one has \E(S)\ = f(\S\). 

Let S a , 1 < a < k arc the C^ 1 = k subsets of N k such that \S a \ = k - 1. Clearly, S 1 U S 2 U 
. . . U S k = N k . Any two components can be potentially related to each other according to (3) only 
if their sets of n indices are related by a permutation. Thus, if k > n the set E(N k ) of independent 
components of A %1%2 "' ln with i a G N k is the union of the k sets E(S a ) of independent components 
with i a e 5" a : E(N k ) = EiS 1 ) U E(S 2 ) U . . . U E(S k ). For the same reason for any two sets Si 
and S 2 one has E(Si) (~l E(S 2 ) = E(Si n S 2 ). Applying the Inclusion-Exclusion Principle [2,7] to 
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the k sets E(S a ) one gets 

\E{N k )\ = \U^E{S a )\= J2 W°)l 

l<a<fe 

- J^ \E(S a )nE(S b )\ 

l<a<b<k 

+ J2 \E(s a )nE(s b )r\E(s c )\- 



l<a<b<c<k 

+ {-l) k - 1 \E{S 1 )n...nE(S k )\, k>n. (14) 

This formula immediately implies 

fc-i 

f(k) = j2(- 1 Y~ 1 f( k -j) c i> k>n - ( 15 ) 

Since this last formula is valid for any k > n one can compute f(k) for k > n from the n values 
/(i) for i = 1, . . . , n. This means that the /(fc) is a polynomial of fc of degree n or less. D 

From the computational point of view it is useful to explicitly compute the above-mentioned 
representation of f(k) through the n values f(i) for i = 1, . . . , n. 

Theorem 2. For an object A lll2 '" ln with arbitrary n and arbitrary compatible linear symmetries 
the number of independent components f(k) for any k > n can be computed as 

n 

f(k) = ^(-l)- 1 /(») CI C n k -_U- (16) 

i=l 

Proof (1). A straightforward way to prove (16) is by induction. For k = n + 1 Eq. (16) coincides 
with (15) and is thus correct. Suppose that Eq. (16) is valid for some k. Let us prove that it is also 
valid for k + 1. Applying Eq. (15) one gets 

k 

/(fc + i) = £(-i) fc -VWCj + i- (17) 

1=1 
Since for all f(i) with n + 1 < i < k Eq. (16) is supposed to be correct, one has 

n k 

f(k + l)=J2(-l) k - i f(i)Cl +1 + Yl (-i) fc_< /WCJ+i 

i— 1 i— n+1 

n /en 



j=n+l 



=^(" 1 )"" i /« 



(-i) fe -Q +1 + J] (-if-^ +1 qc7p. 



j=n+l 



=^Zi-^ l - l f^)Cl +1 ClZl (18) 

i=l 

Here we used the well-known properties of the binomial coefficients [2, 7]. In particular, we used 

k—n 

£(-l) S Cr#i 1+S C:_i +s = l> 0< ? <n<fc. (19) 

s=0 

Therefore, Eq. (16) is correct for any k > n. □ 
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Proof (2). A more elegant proof of Eq. (16) directly follows from Theorem 1: if a function is a 
polynomial of degree n or less, one can take its first n values f(i) for i = 1,2, ... , n and construct 
a polynomial of degree n having these values for i = 1,2, ... ,n (and /(0) = 0) using the Lagrange 
interpolation formula. Indeed, 

(k - 0)(fc - 1) . . . (fc - (i - l))(fc - (i + 1)) . . . (k - n) 
(i-0)(i-l)...(i-(i-l))(i-(i + l))...(i 
which is exactly the Lagrange form of the coefficients of the interpolation polynomial. □ 



(-1)™-' c\ c n k zU = \, n ^, /,••,, h tttt: , ,» . „/ . ( 20 ) 



5 On the Compatibility of the Symmetry Properties 

Up to now we have ignored the question of compatibility of the symmetry properties (3). One can 
easily check if the symmetries are compatible for k = 1. 

Theorem 3. For an object A 1112 '" 1 " with arbitrary n its symmetry properties (3) are incompatible 
for k = 1 if and only if at least one of the two conditions are met: 

(1) for at least one symmetry bi ^ and Q = with 

pi 

c ' = Xl a ^' ( 21 ) 

(2) there exist at least two symmetries for which Ci ^ and c\i ^ 0, and b\jci ^ b\< /cy . 
Proof. For k — 1 any of the symmetries can be written as 

c l A 11 - 1 = b h (22) 

Ci being defined by (21). These properties can be analyzed in the same way as we did for (7). The 
theorem immediately follows from this analysis. □ 

Another two important results are: 

Theorem 4. If symmetries of an object A lll2 '" ln with arbitrary n are compatible for dimension 
k = K , they are also compatible for any dimension k < K . 

Proof. This is obvious since the whole system of symmetry-induced equations for any k < K is a 
part of the corresponding system for k — K. □ 

Theorem 5. If symmetries of an object A lll2 '" ln with arbitrary n are compatible for dimension 
k = n, they are also compatible for any dimension k. 

Proof For k < n this follows from Theorem 4 while for k > n from Theorem 2. □ 

The algorithmic usage of Theorem 5 is clear: one has to check if the symmetries are compatible 
for k = n. 

6 Algorithms to Compute the Independent Components 

Below three different algorithms to find the independent components of an indexed object are 
discussed. 

6.1 Algorithm A 

The first algorithm is a straightforward one. For a fixed k all possible combinations of the numerical 
values of all n indices are substituted into each of the symmetries (3). Thus, one obtains a system 
of linear equations for all k n components of the object. The total number of equations is sk n , 
where s is the number of symmetries (s < n\). Each equation involves at most n\ components. 
Then, the system is solved and the independent components are found explicitly. Counting them 
allows one to get f(k) for that fixed k for which the system was generated. It is clear that if n, k 
or s are large the calculations can be very time-consuming. The only reason to implement such an 
algorithm is the possibility to check better algorithms described below. 
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6.2 Algorithm B 

From the algorithmic point of view, Eq. (16) allows one to calculate f(k) for any k > n as 
soon as one has calculated f(k) for all 1 < k < n. Moreover, having the set of the independent 
components for k — n one can count the number of components among them with indices i a £ N p , 
N p = {1, . . . ,p} for p = 1, ..., n— 1. The number of such components is exactly f(p), and, therefore, 
it is sufficient to have the set of independent components for k = n to compute f(k) for any k. 
Such an algorithm requires solving at most s n n linear equations with at most n\ unknowns in each 
equation and at most n n unknowns in the whole system (note that each of the s < n! symmetries 
(3) generates at most n™ linear equations for at most n\ unknowns). Although, this algorithm is 
better than straightforward calculation of the independent components for some large k > n, it is 
still quite time-consuming for larger n. 

Another point is that this algorithm does not allow to list the independent components and 
to represent the other components as functions of the independent ones for k > n. It is certainly 
possible to augment the algorithm in this direction. However, attempts to do so allowed the author 
to formulate much more efficient algorithm for both computing f(k) and finding the dependencies 
in explicit form. This algorithm described in the next Section. 



6.3 Algorithm C 

Here we suggest a much faster algorithm based on the fact that the sequences of indices of the 
components which could be potentially constrained by a symmetry of the form (3) are related 
with each other by a permutation. This obvious fact has been already mentioned and used above. 
Therefore, we can split the whole set of k n components into such subsets within which the sequences 
of indices are related by a permutation and then generate and solve the corresponding linear 
equations only for the components within each of these subsets. Moreover, one can drastically 
reduce the number of the subsets to be considered since many of them are similar to each other 
(e.g., they can be obtained from each other by changing, say, value 4 for all indices into value 5). 

Let us consider a component A %1%2 '~ ln with some fixed indices 1 < i a < k. Any sequence of 
indices iii2 ■ ■ ■ % n can be characterized by sequence X = (xi, . . . ,£&), where each Xb is the number 
of such i a in iii 2 • • • in that i a — b. Clearly, one has Xi+x 2 + . ■ -+Xk = n with constrains < x\, < n. 
For fixed n and k there are C k+n _ 1 different solutions of this equation with these constrains, and, 
therefore, C^ +n _ l different sequences X. 

It is clear that two components A %ll2 '~ ln and .A-? 1 -? 2 ■■•■'» can be related to each other by a 
symmetry of the form (3), only if both sequences of indices i\i 2 . . . i n and jij 2 . . . j n correspond to 
the one and same sequence X . The number of components A lll2 '" ln corresponding to the same X 
is the multinomial coefficient (n; Xi, x 2 , . . . , Xk) — n\ / (x\\x2\ ■ ■ ■ £fc!) with x\ + X2 + . . . + Xk = n. 

The subsets of the components corresponding to two different X\ and Xi can be treated in 
the same way if X\ and Xi are related to each other by a permutation. Such a permutation 
corresponds just to renaming all indices having, say, value 1 to, say, 5, and so on. The two subsets 
corresponding to two such X\ and Xi have the same number of independent components and the 
same dependence of the other components on the independent ones. Therefore, it is sufficient to 
calculate the dependence of the components only for one X among all of them related to each 
other by a permutation. We will consider only the sorted version Y of X: Y = (yi, . . . ,yk) with 
2/i > 2/2 > ••• > Vk- 

Let Y l be the sequence Y with exactly I nonzero elements: Y l = (j/i, . . . ,yi,0, . . . ,0), \Y l \ = k, 
2/1 > 2/2 > • • • > Vi > 0. It is obvious that 1 < I < min (n, k). Each Y l corresponds to a partition of 
n into I parts. There exists P(n, I) different partitions of this kind. 

Now, let p < I denote the number of distinct values among 2/1, 2/2, • • ■ , yi, and 1 < s m < n, 1 < 
to < p is how many times the value number m (which is one of the p distinct values in Y ) appears 
among j/i, 2/2, • • • , yi- Then the number of different X which can be obtained by permutations from 
Y l isC l k l\/(s 1 ls 2 \...s p l). 
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Combining all the results discussed above one gets the following formula for the number of 
independent components 



f(k) = 



min (n,k) 

E 
1=1 



P(n,l) 



E^ 



v. 



„3\„3\ 



<J | 






(23) 



where the inner sum goes over all the partitions of n into I parts, each such partition corresponds 

s J 2 ,...s J p is calculated for Yj, and g(YJ) is the number of independent 



to Yl, the numbers s{, 
components among n\/(yi\ . 



■ yv- 



components A* 1 * 2 "*" corresponding to Yj: 



0<90l)<^ 



■ yv- 



(24) 



Two identities can be used to check the internal consistency of (23) and (24). First, the total 
number of sequences X for fixed n and k can be calculated from (23) with g(Yj) — 1 an d should 



beCT 



as discussed above. One can sec that this is true: 



min (n,/c) 

E 
i=i 



P(n,l) 

E 



s{\ 









(25) 



Second, for an object without symmetries g(Y 



h =n\ 



n!/(j/i! 



. yi\), and the total number of compo- 



nents for fixed k and n calculated according to (23) should be k n . Indeed, it is also true: 



min (n,k) 

E 
i=i 



P(n,l) 

E 



2/1 1 



■ yv- 



s{\. 



j\ 



c: 



(26) 



Combining (23) and (16) it is clear that in order to calculate f(k) for arbitrary k it is sufficient 
to calculate the number of independent components within the subsets of A* 1 ' 2 '"*™ corresponding 
to Y^i=i P( n i — P{n) different sequences Yj , P(n) being the total number of partitions of n. The 
size of the subsets of A ni2 ' 4 " does not exceed n!. Therefore, in the worst case one should solve 
(n!) 2 linear equations (each symmetry (3) generates at most n! distinct equations, and there are 
at most n! symmetries) with n! unknowns. This is much better that for the algorithms A and B. 

6.4 Reduction of the Number of Linear Equations 

A simple idea allows one to reduce further the number of equations in the system generated by 
the symmetries before solving that system of these equation. One can put the equations into 
a canonical form in which it is trivial to check if any two equations are equivalent (e.g., the 
numerical coefficient at the lexicographically first component should be equal to unity) and retain 
only one among the equivalent equations possibly appearing in the set of generated equations. 
For example, the symmetry T y — T' Jt = for i G {1,2} produces two equations T 12 — T 21 = 
and T 21 — T 12 = which are equivalent and can be considered as one equation. On the other 
hand, the symmetry 3T tJ + AT Jl = 7 for i G {1,2} gives two linearly independent equations 
3T i2 + 4T 2i _ 7 and 3T 2i + 4T i2 _ 7 whose canonica i f orms are different: T 12 + | T 21 = | 

and T 12 + jT 21 = |, respectively. How many equations can be eliminated from the system of 
equations using this simple equivalence test depends on the symmetry properties (note, that even 
two different symmetry properties can produce equivalent equations). This reduction scheme for 
the system of linear equations can be used in all three algorithms described above. 



7 Implementation in Mathematica 

In order to check the performance and cross-check the results all three algorithms A, B, and C 
were implemented in Mathematica with the idea that the best one should be incorporated into the 
package Lino for calculations with indexed objects [9,8]. The main parts of the implementation 
are: 
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(1) two routines Def Object and Def Symmetries allowing one to define objects with arbitrary 
symmetries, 

(2) routine ConstrainComponents which explicitly generates and solves the linear equations for 
individual components induced by the symmetries for some fixed k, 

(3) routine GuessPolynomial implementing algorithm A by calling ConstrainComponents for a 
sufficient number of different values of k to check if a polynomial of degree n or less can be 
fitted to the results, 

(4) routine CountlndependentComponents implementing algorithm B by calling ConstrainComponents 
for k = n and analyzing the resulted independent components to compute f(k) from (16), and 

(5) routine ListlndependentComponents implementing algorithm C and providing for any k both 
f(k) and, if requested, a list for the independent components and the dependence of the other 
components. 

All the routines allow one to control all the steps of the corresponding algorithms and, if desired, 
provide the user with various additional information. The implementation consists of about 1000 
lines of Mathematica top level code and is available form the author upon request. 

To give a practical example let us consider the covariant Riemann tensor Rijki with n = 4 and 
with its four symmetry properties 

ti-ijkl — ftklij > 

ti-ijkl ttjiklj 

K-ijkl ttijlkj 

Rijkl + Riljk + Riklj = 0. (27) 

The well-known result [11] for the covariant Riemann tensor is f(k) = j^k 2 (k 2 — 1). 

Algorithm A (GuessPolynomial) explicitly computes the independent components subsequently 
for k = 1, k = 2, . . . , k = 6 to verify that a polynomial of degree 4 or less really fits the results. 
For example, for k = 6 this requires solving of a system of 4 x 6 4 = 5184 linear equations (2526 
distinct ones) with for 1296 unknowns. This is quite a heavy task even for a modern PC (note that 
the system is undcrdctcrminatcd and should be solved exactly). 

Algorithm B (CountlndependentComponents) requires solving a total of 4 x 4 4 — 1024 linear 
equations (504 distinct ones) for 4 4 = 256 unknowns among which 112 components turn out to 
be zero and another 124 turn out to be functions of 20 independent components: -R2121J -R3121, 

-R3131, -^3221, -R3231! ^3232, ^4121, -R4131j -^4141, R4221, -^4231; -^4232, #4241, #4242, #4231, #4331, 

#4332, #4341, #4342, #4343- This list allows one to conclude that /(l) = 0, /(2) = 1, /(3) = 6, and 
/(4) = 20. These set of values allows one to compute the above-mentioned result for f{k) directly 
from (16). 

Algorithm C (ListlndependentComponents) required solving of #(4) = 5 systems of linear 
equations: 

(1) a system of 4 equations (only one distinct equation) with 1 unknown (corresponding to Yf = 

(4,-..)), 

(2) one system of 16 equations (9 distinct ones) with 4 unknowns (corresponding to Y 2 — (3, 1, ...)), 

(3) a system of 24 equations (10 distinct ones) with 6 unknowns (corresponding to Y 2 = (2, 2, . . .)), 

(4) a system of 48 equations (24 distinct ones) with 12 unknowns (corresponding to Yj 3 = (2, 1,1,.. .)), 

(5) a system of 96 equations (44 distinct ones) with 24 unknowns (corresponding to Yj 4 = (1,1,1,1,...)). 

Clearly, algorithm C produces the same results for f(k) as the other two algorithms, but requires 
much less resources. This demonstrates the efficiency of algorithm C. It is planned to include 
algorithm C into the next release of Lino. 

8 Concluding Remarks 

Certainly, algorithm C can be further improved in certain cases if the structure of particular 
symmetries are taken into account. Up to now the algorithms does not account for any properties 
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which the symmetries may have. Here one can use the group-theoretic approach for manipulations 
with indexed objects developed in [14,5,6, 13]. Certain further improvement could be achieved if 
an algorithm to generate only distinct equations could be found. However, it is doubtful that such 
an algorithm would be computationally cheaper as the currently used algorithm to find and to 
drop equivalent equations before solving the system of linear equations (sec, Section 6.4). It is 
interesting also to check if the results of Section 5 can be improved so that the incompatibility of 
symmetries (3) could be seen in an easier way. Another interesting question is whether, for arbitrary 
symmetries, one can express algorithmically the combinatorial "finger exercises" allowing one to 
derive /(2) = 1, /(3) = 6 and /(4) = 20 for the covariant Ricmann tensor as given, for example, 
in Section 92 of [11]. 

The considered form (3) of the symmetries does not allow us to consider some important cases. 
For example, the definition of a symmetric trace-free (STF) tensor required that a contraction of 
■^-i\ii...i l with the Kronecker symbol Si a i b vanishes for any a and b. Such a symmetry cannot be 
written in the form (3) and is out of the scope of this paper. On the other hand, STF tensors plays 
very important role in modern physics [1,4, 10] and it is important to have efficient algorithms to 
store them and manipulate with them. It can be demonstrated that the main results of this paper 
can be also used for symmetries involving contractions with objects each component of which has 
some numerical value. That is, one can consider symmetries of the form 



p 

E 



a, J2 B ili2 ... im AT* =6, (28) 



where clj and b are numbers, i\, ii, ..., i m are dummy indices over which the contraction is 
performed, -B^^...^ is a number for any values of its indices (this can be, e.g. the Kronecker 
Sij or the fully antisymmetric Levi-Civita symbol Eijk, or anything else), and Jfj is an arbitrary 
permutation of n indices containing m dummy indices i\, %2, . . . , i m an d n — m free ones. This 
case will be treated in a separate publication. 
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